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ABSTRACT 



We analyse subarcsecond resolution interferometric CO line data for twelve sub-millimetre- 
luminous (Ssso^m > 5 mjy) galaxies with redshifts between 1 and 3, presenting new data for 
four of them. Morphologically and kinematically most of the twelve systems appear to be major 
mergers. Five of them are well-resolved binary systems, and seven are compact or poorly resolved. 
Of the four binary systems for which mass measurements for both separate components can be 
made, all have mass ratios of 1:3 or closer. Furthermore, comparison of the ratio of compact to 
binary systems with that observed in local ULIRGs indicates that at least a significant fraction of 
the compact SMGs must also be late-stage mergers. In addition, the dynamical and gas masses 
we derive are most consistent with the lower end of the range of stellar masses published for 
these systems, favouring cosmological models in which SMGs result from mergers. These results 
all point to the same conclusion, that likely most of the bright SMGs with Ljr > 5xl0 12 L Q are 
major mergers. 



Introduction 



Surveys f or dust-obscured g alaxies at high redshifts carried out with SCUBA ((Holland et al.l Il999l ) 



(5Wi / /rn >5mJ 
booa IPope et al 
(~1O 13 L0), intensely star- forming ( 



20051: 



Weiss et al 



Greve et al. 


2004: 


Coppin et al. 


2005; 


Chapman et al. 



20091 ) . This population c onsists of dusty, gas-rich, very luminou s 
^Moyr -1 ) galaxies (e.g. iGreve et al.lbooalchapman et al.ll2005l) . 



The importance of these SMGs is underlined by studies of the extragalactic far-IR/sub-mm background, 
which shows that about half of the c osmic energy density originates in distant, dusty starb ursts and 
AGN (buget et al.lll99d: IPei et al.l |l999|) with ULIRGs contrib uting -20% jMagnelli et alJboiol). Near-IR 
( Swinbank et al. l2004h . optical and X-ray ([Alexander et al . 2005), and mid-IR ( Lutz et al . 2005; IValiante et al 
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2007 ; Pope et al. 2008 ; Menendez-Delmestre et al. 20091) spectroscopic work has shown that starbursts are 
the primary power source, with AGN present in a large fraction but not making a dominant energetic 
contribution. 

A key question is what triggers the extraordinary luminosities of these galaxies. A number of authors 
have suggested that at least some SMGs may be interacting or merging systems, based primarily on irregular 



or complex morpholo gies seen in rest-frame UV HST images (jChapman et al.ll2003l l2004HSmail et al.ll2004 ; 
Swinbank et al. [20061) . However, two caveats need to be kept in mind regarding these results. UV emission 



is strongly affected by extinction, and hence UV images may not reliably track the u nderlying galaxy. This 
is demonstrated by the radio observations carried out by Chapman et al. ( 2004 ) and Biggs fc Ivison ( 20081 ) 



(a key result of these studies was that the star formation in SMGs is more extended than in local ULIRGs, 
an issue which we discuss in more detail later) , which show morphological differences between radio and UV 



Lotz et al. 


2006; 


Conselice et al. 


2008) 



the consequen ce of clumpy star formation in a kinematically more coherent structure such a s an underlying 



rotating disk (iForster-S chrcibc r et al 



Ivison et al 



(200 



2006 
I 



jool) were able to show that SMGs have an increased likelihood 



2009; iGenzel et all 120061; Elmegreen et all 120071) . or strong and 



structured extinction. _ 

of having close neighbours, suggestive of physical associations. Swinbank et al. ( 20061 ) obtained integral field 
near-infrared spectroscopy data of 6 SMGs, which showed a significant fraction to consist of kinematically 
distinct sub-components, interpre ted as an indicat i on of merging. However, as with the HST data discussed 
above and as also pointed out bv ISwinbank et al.l (|2010l ). these rest-frame optical and UV data suffer from 
extinction limitations and potential biases due to outflows. Ideally we need spatially resolved maps and 
kinematics of a tracer that is not affected by dust exti nction to properly assess m orphologies and kinematics. 
Using high resolution millimetre CO line observations. iTacconi et al. | |2006ll2008l ) find SMGs to display either 
double/multiple morphologies with complex gas motions indicative of on-going mergers, or to be extremely 
dense, compact, rotating disks. The high gas and matter densities of the latter argue for them to be end-stage 
major mergers. However the spatial resolution achievable with current facilities is insufficient to definitively 
distinguish between coalesced major mergers and isolated disk galaxies. At present it is thus still unclear 
whether SMGs form a homogeneous class of galaxies, or whether their high luminosities can be achieved 
through a variety of pathways. Illustrative of this on-going discussion are the deb ates about the nature of 
two of the first SMGs to be discovered: SMM J02399-136 and SMM J 140 11 +0252 dlvison et al.lll998l l2000[ 

IT" 



Smail et al 



2002; 



Genzel et al 



2003 



Ivison et al 



2010; 



Downes fc Solomon 



2003 



Smail et al 



200 



5). 



The aim of this paper is to critically examine the hypothesis that SMGs are mergers. We use observations 
of the CO line, which is unaffected by extinction or outflows, to trace the molecular gas, and combine 
both morphological and kinematic aspects in the analysis. We present new data for four SMGs. One of 
our sources has not been observed in CO before (SMM J105141), for the three other systems these new 
observations provide data of significantly increased resolution. This brings the number of SMGs observed 
to sub-arcsecond resolution from 8 to 12 and thus provides a significant increase in the total sample of 
high resolution CO SMG observations. We complement these new data with previously published high 
resolution CO data for eight other SMG s jGreve et al.ll2005t Elainline et alJkood : ITacconi et al.ll2006l l2008t 



Ivison et al 



201 



Tacconi et al 



0; 



(200 



■Bothwell et al.ll201(i| ). resulting in a much- improved statistical basis compared to e.g. 
^" 2008 ). Our work is distinctive to these preceding works, the aim of which was a general 



characterisation of the SMG population, not only in regard to sample size, but also objective. Here we tackle 
a specific, heretofore unsolved key issue - namely, whether or not SMGs are major mergers. 



In Sj2] we detail our observations and data reduction. In Sj3] we explain the analyses employed, and 
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present results for individual systems. We then investigate merger mass ratios, and sizes of our sample 
galaxies. In ij?]we test stellar mass estimates from the literature, and discuss the implication of our results 
for cosmological models. We then bring these results together to discuss the trigger mechanism for the 
exceptional luminosities of SMGs, before concluding in 33 



2. Observations & Data Reduction 



The observations of four new SMGs were carried out in wint er 2007/2008 with the IRAM Plateau de 
Bure interferometer, which consists of six 15m-diameter telescopes ([Guilloteau et al. Bl992h . We observed the 
four sources, SMM J09431, SMM J105141, SMM J131201, and HDF 169 under excellent weather conditions in 
the extended A-configuration, which has baselines up to 760 meters. The resulting spatial resolutions ranged 
from 0.4"-0.6"(see Figs. HE]). For the two z~3 sources, SMMJ09431 and SMMJ131201, we observed the 
CO J = 6-5 line, which is redshifted into the 2mm band. For the two sources at z — 1.2, SMMJ105141 and 
HDF 169, we observed the CO J = 4-3 line in the 1mm band. The observations were all made with the 
new generation, dual polarisation receivers, with system temperatures of ~150K and — 200K (SSB) in the 
2 mm and 1 mm bands, respectively. The full available correlator capacity wa s configured to cover 1 GHz in 
each polarisation, which adequately covers the broad CO lines seen in SMGs |Genzel et al.ll2003t Neri et al 



20031 ; iGreve et al.l 120051) . The total on-source integration times were 9.3 hours for SMMJ09431, 12.7hours 
for SMM J105141, 14.4hours for HDF 169, and 20.1 hours for SMM J131204. 

We calibrated the datasets using the CLIC program in the IRAM GILDAS package. Passband cali- 
bration used one or more bright quasars. Every 20 minutes we alternated source observations with a bright 
quasar calibrator within 15 degrees of the source. Our absolute positional accuracies are ±0.2" or better. 
The absolute flux scale was calibrated on MWC 349 (1.35 Jy at -156 GHz and 1.6 Jy at -208 GHz). After 
flagging bad and high phase noise data, we created data cubes in the MAPPING environment of GILDAS. 
Final maps were cleaned with the CLARK version of CLEAN implemented in GILDAS. The absolute flux 
scale is better than ±20%. 



3. Results 

Integrated spectra, integrated line flux maps, position-velocity diagrams, and velocity & velocity dis- 
persion maps of the individual systems are displayed in Figures [TJ [2j [3l & [4] Measured and derived source 
properties are listed in Tables Q] [3] & [2j 



3.1. Analysis 



Our analyses generally follow th e procedures of Tacconi et al. ( 2006 . 20081) ; for the computation of cold 
gas masses, we use the approach of Tacconi et al. ( 2010l ). Flux maps were produced by integrating the 
line emission, determining the line extent via visual inspection of spatially-integrated spectra. Fluxes are 
extracted from the line cubes by integrating over the area with significant emission. Given that the compact 
galaxies comprising a significant fraction of our sample here are in most cases only marginally resolved 
with the ~ 0.5-1" spatial resolution available, we determine source sizes by fitting circular Gaussians to the 
systems in the ww-plane using the task UVFIT in the IRAM data reduction package GILDAS, and taking 



-4- 



the Gaussian FWHM utI to be the best estimate of the source's diameter. In the two individual galaxies of 
the four new systems where a Gaussian fit was not possible (most likely due to insufficient signal-to- noise) , 
the beam FWHM was used as an upper limit to the source diameter. We note that these size measurements 
may slightly underestimate the actual size of the galaxies, since there may be faint extended flux below our 
sensitivity limit. We do not assign a formal uncertainty to this, since, if such a putative faint emission region 
exists, we do not know anything about its size or morphology. However the uncertainties due to this are 
unlikely to be significant. 

Where possible, velocity and velocity dispersion maps were derived by fitting Gaussians to the line 
profiles. If the spectral signal-to-noise did not permit this, we produced moment maps of different velocity 
ranges, which were then overlaid as RGB-images. 

Since the dispersions of our galaxies are comparable to or larger tha n the rotational velocities, dynam- 
ical masses have been derive d using the 'isotropic virial estimator' (e.g. Spitzer 1987 ; Tacconi et al. 20081 : 
Forster-Schreiber et al.ll2009l) given in M & by 



dyn,vir 



= 2. 



x W 5 Av 2 FWHM R 1/2 



(1) 



where Avfwhm (kms -1 ) is the integrated line FWHM, and R\/2 (kpc) is the the half-light radius (the 
radius within which half the total flux is emitted, here assumed to be equivalent to the Gaussian HWHM 
measured in the iw-plane). We note that the scaling facto r appropriate for a rotating disk at an average 



inclination is a factor of ~1.5 smaller (Bothwell et al 



2010h . i.e. if the ab ove formula were applie d to a disk 



galaxy, it would somewhat overestimate the dynamical mass. However, ICappellari et al.l fj2009f ) find that 
this estimator agrees well with masses derived from more detailed Jeans modelling for a sample of massive 
z ~ 2 early-type galaxies. 

In order to derive molecular gas masses, we convert th e measured line fluxes to CO(l-O) fluxes (in 
Jykms" 1 ) using the results of SMG SED modelling (Fig. 3 of IWeissl[2007l yielding £?co(6-5)/£co(i-o)=13, 
Snnu-z\ /noh-()\ = l O, S rmt*--?.) I Sr.q^ -ni = 7, Sco(2-i)/<Sco(i-o)=3; the same numbers were also used by 



Bothwell et ail l2010l and lCasev et alJI'2009bl ). We assign an uncertain ty of 30% to these conversions based 



on the differences between the individual SED models of IWeissI (|2007l) : in the four cases where we have flux 
measurements from two different transitions (Table [l} we find the flux ratios to be consistent with these 
conversion factors. We do not calculate gas masses based on CO (7-6) fluxes, since the conversion factor is 
u ncertain for this transi tion. We then calculate the CO(l-O) line luminosities (in Kkms -1 pc 2 ) using Eqn. 3 
of ISolomon et al.l |l997h : 



CO(l-O) 



3.25 x 10 7 S CO(1 _ 0) A^- 2 (1 _ 0) rest D\{\ 



(2) 



with Sco^V m Jykms , vco(i-o),rest m GHz, and Dl in Mpc. The CO(l-O) luminosities are then used 
to estimate the total gas mass via their Eqn. 19: 



Ms, = ah'. 



CO(l-O) 



(3) 



using a = 1 / K km s 1 pc? as measured for ULIRGs ( Downes fc Solomon 19981 for a disc ussion see also 
Genzel et allboioh . For an in-depth discussion of this we refer the reader to Appendix A of iTacconi et al.l 
(200 81); we would like to stress that our choice of a is motivated by the surface densities of SMGs (cf. Table 



10 of Tacconi et al. 20081 ) rather than any a priori assumption about a similarity between local ULIRGs and 
SMGs. Our choice of conversion factor is here further supported by the fact, that if for our SMGs the Milky 
Way conversion factor was adopted, in all but one case the gas mass would, rather implausibly, be equal 
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to or larger than the dynamical mass. Finally, we apply a factor 1.36 correction to account for interstellar 
helium. We note that since the fluxes were measured in fixed apertures, and also some flux may be resolved 
out or be below our detection limit, the gas mass is likely to be somewhat larger than our estimates. 



3.2. Results For Individual Objects 

The detailed analyses of the molecular CO line emission in this Section reveal morphological and kine- 
matic evidence that at least eight of the twelve SMGs studied in this paper are on-going mergers, and the 
remaining four are either disk galaxies or coalesced, late-stage mergers. 



3.2.1. New Observations 



SMM J09431 (z = 3.35) - a review of previous observations is given by Tacconi et al. ( 2006b . who them- 
selves presented CO(4-3) data. The system consists of two galaxies, H6 and H7, separated by ^30kpc in 
projection, thus physically related and perhaps in an early stage of merging. Our new CO (6-5) data show 
a clear line spectrum for the so uth-west galaxy H7, with a gl obal velocity gradient from -50 to +70kms -1 , 
probably indicative of rotation. Forster-Schreiber et al. ( 20091) find that the commonly adopted boundary be- 
tween 'dispersion-dominated' and 'rotation-dominated' systems of v ro t ~ &o (where v ro t is the circular/orbital 
velocity and o~$ is the intrinsic velocity dispersion) translates into 0~i nt /v o b s ~ 1.25 (where v b s is the observed 
line-of-sight rotation and <r mf is the integrated line FWHM divided by 2.355). H7 has <Jmt/v bs > 1-5 nearly 
everywhere, i.e. appears dispersion-dominated. For the north-east galaxy H6, we detect a 4er line shifted 
about +210kms~ 1 relative to H7. Its width is 180 ± 65kms _1 . There may also be broader and more 
redshifted emission extended over a larger region around H6 and towards H7, but at low significance. When 
deriving gas masses, the measured fluxes an d sizes are corrected by a factor 1.3 for gravitational lensing 
effects, using the model of Cowie et al. ( 20021 ). 



SMM J131201 (z = 3.41) - This system shows complex spatial and velocity distributions. There is no 
evidence of ordered rotation, although the extensions to the east and north-east are redshifted with respect 
to the emission peak. The substructures, which may be clumps, outflows, or tidal tails, and the overall lack 
of rotatio n most likely p oint towards an advanced, pre-coalescence merger stage. This conclusion is also 
drawn bv lHainline! (j2008h . 



SMM J105141 (z = 1.21) - There are two galaxies separated in velocity by ~400kms _1 and spatially by 
^5 kpc, with no indications of disturbances. The south-west galaxy is about a factor 3 less luminous than 
its companion, but nearly equal in dynamical mass. The line profiles are well-resolved. The velocity field of 
the north-east system is indicative of ordered rotation, with velocities ranging from 60 to 180 km s" 1 . A earn 
we find that generally 0~i n t/v o b s > 1.5. It is likely that we are seeing two galaxies engaged in the early stages 
of merging, appearing close only in projection. 

HDF169 (SMM J123634+6212, z = 1.22) - The integrated spectrum of this compact galaxy is double- 
peaked, and its P-V diagram resembles a rota tion curve, both indic ative of a compact rotating system, 
has a mass surface density of ~ 10 3 3 Mq pc~ 2 . 



It 



Tacconi et al. ( 2008h argued that the very high volume and 



surface densities characterising compact systems such as this could only plausibly be achieved during the 
final stages of a major merger. However we feel that this argument alone here is insufficient to judge with 
any degree of certainty whether this system is a quiescent disk galaxy or a late-stage merger. 
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3.2.2. Re-Analysis Of Existing Observations, & Published Results 



W e have re-analysed previous observations of HDF 242 (SMM J12370 7+614, z = 2.49. ITacconi et al. 
2006) and N2850.4 (SMM J163650+4057 and SMM J16368+4057, z = 2.38. Eacconi et al.ll2006l) using the 



methods described in Section [3~T1 (see Table H). W e al so complement our da ta with published results from 



CO line observations bv ITacconi et all (|2006l . 120081 ) and lBothwell et al.1 (j2010h for five other systems; HDF 76 
(SMM J123549+6215, z = 2.20), N2850.2 (SMM J163658+4105, SMM J16370+4105, and SMM J16366+4105, 
z = 2.45), SMM J044307+0210 (SMM J04431+0210, z = 2.51), HDF 132 (SMM J123618. 33+621550. 5, z + 2.00), 
and Lockman38 (SMM J105307+572431.4, z — 1.52), where possible re-calculating derived quantities to bring 
them into a consistent methodolog ical framework. F inally, we include results from new multi-wavelength 
observations of SMM J02399-0136 (jlvison et al.ll2O10l) . The published results on the kinematics and mor- 



phologies of these systems are summarised below: 

SMM J044307+0210 - This lensed galaxy appears compact and has a double-peaked line profile with a 
FWHM of - 350 km s" 1 . 

HDF 242 - a wide binary with a separation of ~22kpc. 

N2850.4 - The CO(3-2) data show a double-peaked line profile and a separation of ^3.3kpc between 
the peaks of the red and blue channel maps. The integrated emission in the higher resolution CO(7-6) data 
directly shows two peaks separated by ~3.7kpc. This system is therefore interpreted as a close binary. We 
infer a mass ratio of 0.68 from the flux ratios of the two components. 

HDF 76, N2850.2, & HDF 132 - All compact sources with double-peaked line profiles. They could be 
either coalesced major mergers, or massive di sk galaxies. They have extremely high m ass surface densities 



of 10 3 - 5 , 10 4 - 5 , and lO 4 O M pc- 2 respectively (|Tacconi et al.ll2008l : iBothwell et al.ll2010D. at least an order of 



magnitude larger than those m easured in the UV/optically-selected SINS galaxies (jForster-Schreiber et al 



2006 



20091 ITacconi et al.ll2008f) and comparable to those measured in the central nuclear gas disk of Arp 220. 



Tacconi et al 



(2008) argue, that to achieve this, the torques provided by a highly dissipative, major merger 



are necessary. 



Lockman 38 - IBothwell et al.1 (|2010l ) interpreted this object as a very late-stage merger, having two 
kinematically distinct components that are not spatially resolved. 



SMM J02399-0136 - New multi- wavelength observations (jlvison et all l2010h reveal that this system 
consists of at least two major components (most likely three, two of which host a BH); it is therefore 
interpreted as an on-going merger. 



3.3. Mass Ratios 



Ivison et al 



Ivison et al.l (20021) first noted a preponderance of radio doubles among SMGs. With a larger sample, 



(200 



3) were able to demonstrate that the likelihood of SMGs having companions within a few 



arcseconds is significantly larger than would be expected by chance. This overdensity suggests that at least 
a significant fraction of SMGs may currently be undergoing mergers. We are able to refine this picture 
substantially by investigating the mass ratios of the merging components. Five of the twelve SMGs in our 
sample consist of spatially separated interacting galaxies, and we are able to measure mass ratios for four 
of them (Table [2 SMMJ09431, SMMJ105141, HDF 242, N2850.4; no mass estimates for the individual 
components of SMMJ02399 are published). In all cases the mass ratios are consistent with being major 
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(mass ratios of 1:3 or closer to unity) mergers. It must be kept in mind that if these are mergers between 
disk galaxies, the dynamical masses may be over- or underestimated due to the dependence on viewing angle, 
and hence the uncertainties on these mass ratios are somewhat larger than the propagated errors stated in 
Table[5J All measured spatial separations (<30kpc) are well within the range where they would be deemed 
gravitationally int eracting, and the masses of ind ividual components are >1O 1O M0, much larger than clumps 
in a single galaxy (jForster-Schreiber et al.ll2009l ). 

Bot hwell et al. I d2010l) outline, one can also make a strong case for Lock- 



We furthermore note that, as 



man 38 and HDF 132 being mergers; they show clear indication of two kin ematically distinct components, 
which are however not individually spatially resolved. Bothwell et al. ( 2O10l) calculate the mass ratios of the 
two subcomponents to be 0.24±0.18 and 0.68±0.14, respectively. However, since these two systems cannot 
unambiguously be classified as mergers, we do not include them in Table [2] 



4. Interpretation and Discussion 
4.1. Sizes 

Size measurements of SMGs based on CO interferometric and mm-continuum data have indicated that 
these galaxies are comparatively comp act, with median diameter <4kpc (FWHM of Gaussian fitted to 
source in wu-plane; Tacconi et al. 20061) . O n the other hand , based on the diameter at the 3 a detection level 
of MERLIN and VLA radio observations, Chapman et al. ( 2004) find the st ar formation to be 'occurring 
on ~10kpc scales'. However. iTacconi et al.l (|2006l ) and iBiggs fc Ivisonl (|20081 ) point out that the size mea- 
surements from the MERLIN and VLA data are consistent with the CO size when a comparable method 
(deconvolved major axes of Gaussians fitted to the emission) is applied. They find a median diameter of 
•~5kpc for their sample. We can compare radio and C O size measurements individually for two sources 
in our sample. For SMM J105141 IBiggs fc Ivisonl (|2008l ) measure an overall radio size of 0.61"x0.44", or 
5.1 kpc x 3.7kpc; we measure a CO(4-3) FWHM of 3.5kpc for the main component of this two-galaxy sys- 
tem. For HDF 169, we find a CO(4-3) FWHM of 4.3 kpc and measure 4.5 kpc from an unpublished MERLIN 
radio map. Thus there app ears to be good agreement between CO(4-3) and radio extents in these two cases. 
Younger et al. I (120081 120101) measure the scale of far-infrared emission of four SMGs to be ^5-8 kp c. This 
is somewhat larger than what we find for the sample discu ssed here (although Younger et al.ll2008l say that 



their results are 'qualitatively consistent with the results of lTacconi et al.ll2006F ). The most plausible expla- 
nation for the difference - besi de the possibility that CO li nes and far-IR continuum trace slightly different 
things - lies in the fact that the lYounger et al] (|2008l l2010h SMGs probe the brightest end of the luminosity 
distribution, with SWo/zm ^ 15 mjy. 

For our augmented sample, which contains all SMGs for which sub-arcsecond resolution interferometric 
CO observations are available (which are not strongly lensed), we have measured FWHM sizes for 11 galaxies 
(measuring each galaxy individually in the case of wide binaries). We find average a nd median FWHMs o f 
4.5 and 4.3 kpc, with a 1-tr distribution of 2.8kpc, slightly larger than the results of lTacconi et al.l (12 006). 
We a re unable to identify a plausible explanation for this, other than small number statistics ( Tacconi et al.l 
I2OO6I measure sizes for five systems, and upper limits for a further four; here we have size measurements for 
twelve systems). We also obtain upper limits of 3.4, 6.5, and 2.8 kpc for three other sources (Table [1]). 

We note that here we are collating measurements fr om a range of CO tr ansitions, from (7-6) to (1-0), 
which may trace physically different environments. As Tacconi et al. ( 20081 ) elaborate, transitions up to 
J « 5-6 should all provide reasonable tracers of the overall molecular gas. However, the extent to which, say, 
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C0(2-l) and CO(6-5) trace the same gas, is still on-going (see also Weiss! 2007 ; Narayanan et al. 2010b ). 
That the CO(7-6) may trace denser gas might explain the comparatively small source sizes for N2850.2 and 
N2850.4, and hence be at least pa rtially responsible for the small dynamical masses we derive for these two 
galaxies (see also the discussion in lDanielson et al.ll2010l ). In the one case where we have a size measurement 
from two different transitions (SMMJ09431 SW; J = 4 & J = 6), the two numbers are consistent within the 
uncertainties. Restricting to size measurements from J < 6 observations yields average and median FWHMs 
of 5.2 and 5.0 kpc, with a 1-a distribution of 2.6 kpc. 



Thus both radio and CO measurements are consistent with SM Gs having FWHM sizes of ^5 kpc. 



4.2. Stellar Mass Measurements, and their Implications for Cosmological Models &: the 

Co-Evolution of Black Holes and their Hosts 

In this Section, we use our dynamical and gas masses to test stellar mass estimates. This is motivated 
by the large systematic uncertainties associated with stellar mass measurements, and the implications that 
a reliable determination of the stellar mass has both for a critical assessment of competin g cosmological 



mode ls of SMGs, as well as for the co-evolution of central black holes and their hosts (e.g. lHainline et al 



20101)- 



4-2.1. Stellar Mass Measurements 



Hainline ( 20081) (see also Hainline et al. 2010 ) derived the stellar masses of a large sample of SMGs, 



employing Spitze r obs ervations to cov er a significantly larger s pectral range than earlie r attempts (e.g. 
Smail et al.ll2004l ). The Irlainlinel (|2008[ ) data set was also used bv iMichalowski et al.1 |2009l ) to deri ye stellar 
masse s. But despite using an identical observational data ba se, their results ar e very different: lHainline 
( 2008 ) find a median stellar mass of 6.3-6.9 x 10 10 M Q , whereas Michalowski et al. ( 20091) arrive at a median 
of 3.5 x 10 u Mq. Therefore, testing these results is essential and the next key step towards deriving the 
stellar mass content of SMGs reliably. 

Table [3] brings our d ynamical and gas masses together with the stellar masses from lHainline! ( 2008 ) and 
Michalowski et al. ( 2009! ) for t he seven sys t ems i n common, and includes gas masses for two more galaxies 
based on flux measurements of Greve et al. ( 20051 ). We also show the stellar masses derived by Tacconi et al. 
(|2008l ). As noted in S j3.1i our dynamical masses may be underestimated due to undetected faint extended 
flux leading to somewhat smaller HWHM UV . However, since this effect would be small and also the gas 
masses correspondingly larger, this is very unlikely to qualitatively impact on our conclusions here. We 
furthermore note that we have neglected b oth dust and dark m atter in our mass budget. Dust is unlikely 
to contribute more than a few percent, and lGenzel et al.l (|2008f ) find the dark matter content in the central 
few kpc of high-z galaxies to be 10-20%. Thus, whilst not negligible, these will not have a significant impact 
on the mass budget, and we h ave opted not t o incl ude it her e for sa ke of simplicity. We have converted the 
stellar mass est imates from a lSahoeterl (Il955l) to a | Chabrierl ([2003) IMF by reducing them by a factor 1.7 



(jHainlind 12008). This is suppor ted by 



Bastian et al.l (120101) . who find no conclusive evidence that an IMF 
different to the locally a pplicable Chabrierl J2003l) orlKroupal ( 2001 ) forms should be in place at high redshifts. 
Cappellari et al. (2009) also find that a lchabrierl (j2003h instead of a Salpeter ( 1955 ) IMF is required to bring 



the stellar masses derived for a sample of z ~ 2 early-type galaxies via SED modelling in agreement with 
their dynamical masses calculated from Jeans models. 
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We find that (M gas + M±)/Mdy n for the seven g alaxies fo r which all three mass estimates are available, 
has an average (median) of 0.8 ( 0.9) for the adjuste d Hainlinej ( 20081 ) stellar masses (with a 1-er distribution 
of 0.5); whereas for the adjusted iMichalowskiet al. ( 2009) m asses, the average (median) is 2.1 (2.1) , with a 



1-cr distribution of 1.2. We conclude that the 
IMF, are consistent with our independent mass measurements (Fig. [6]). 



HainlineT i 20081) stellar masses, adjusted to a IChabrierl d2003l) 



4-2.2. Implications for Cosmological Models 



Two different theoretical approaches to reproduce the sub- millimetre galaxy population in simulations 
are currently being pursued. In semi-analytic models (e.g. Swinbank et al. 2008h . the high sub- millimetre 
luminosities are achieved throu gh bursts of star fo rmation induced in gas-rich major mergers (although 
disk-instabilities may also occur, Bower et al. 20061) : whereas numerical simulations have also invoked 'cold 
accretion flows'. In the latter SMGs are assumed to be massive galaxies sitting at the centres of large 
potential well s, constantly form ing stars at high rates sustained by smooth infall and accretion of gas-rich 
satellites fe.g. lDave et al.ll2010l ). 



Both models have shortc omings (such as cold accretion flows being unable to explain the short SMG 
lifetimes, iTacconi et al.l 120081 ) and it is likely that both m odels are fundamentally limited by our current 
understanding of the role of gas physics inside the haloes. As lBenson fc Bowerd (|2010f ) show, implementation 
of cold mode accretion into semi-analytic models does not translate into a significant difference in cosmic 
star formation history compared to previous models; this may indicate that a more crucial aspect is how the 
stability of the cold gas inside the haloes is modelled. Equally, the assumptions underlying the numerical 
simulations have not been fully tested as it is computationally challenging to run large numbers of galaxies 
in these models to z = to compare their properties to the observations of local galaxy populations (e.g. 
luminosity function, colour-luminosity relation) to confirm that they reproduce these. Clearly, a detailed 
discussion of this is beyond the scope of this paper. Instead, here we will look at the stellar masses pr e dicted 
by the current implementations of each approach, and how these compare to the adjusted Hainlinel j2008i ) 



masses. 

Currently, both modelling approaches have the difficulty that if the number dens ity of SMGs is ma tched, 
then with a standard IMF, they underpredict t he luminosities /SFRs and vice versa (jPave et al.lboiol) . This 
may be due to systemati cs in the simulations ([Swinbank et al.1 120081 ) or the inferred properties of observed 
SMGs (jDave et alJ 120101 ) . A key difference between the models is the stellar mass required to match the 
observed number density. In the cold accre tion flow models , SMGs are the most massive galaxies and 
have a median stellar mass of 2.7 x lO n M0 (jDave et al.l 120101 ) . In the semi-analytic models, the declining 
stellar mass function, major merger rate and brief SMG phase pu sh the predicted stella r mass down and 
the median is an order of mag nitude lower at ~2.lt? n x 10 10 M ? (ISwinbank et al.ll2008l) . In Sectio n J4~2~T1 
we showed it is the lower stellar masses, calculated by Hainlinel ( 2008 ) and corrected to a Chabrier (2003) 



IMF, that are physically more plausible for SMGs. The adjusted median stellar mass for their 64 SMGs is 
3.9 x 1O 1O M0, consistent with what is required by semi-analytic models. We therefore conclude that, although 
the associated uncertainties are certainly non-negligible, SMG stellar mass measurements strongly disfavour 
cold accretion flow models, and are compatible with hierarchical merging models. 
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4.3. What Triggers The Exceptional Luminosity of SMGs? 



It has previously been conjectured, based primarily on radio and rest-frame UV morphologies, that the 
high luminosities of SMGs are triggered by merger activity. Our analysis in Section [3] of the CO emission 
shows convincing evidence that at least a majority of our twelve SMGs are indeed mergers. In particular, we 
have found that all binary SMGs where we could derive mass ratios the merger partners are near-equal in 
mass, i.e. they are major mergers. While it is thus beyond doubt that some SMGs are triggered by major 
mergers, the question of the nature of the compact, non-binary SMGs comprising a significant fraction of 
the overall SMG sample is yet to be settled conclusively - are they end-stage major mergers, or isolated disk 
galaxies? 

In the following we will develop the argument that they mostly are end-stage major mergers, since a 
major merger will be detectable as an SMG for a significant amount of time spent in the 'final coalescence' 
phase - and thus that the large fraction of observed binary SMGs imply a correspondingly large number of 
single SMGs which are coalesced major mergers. 

An illust rative compar i son ca n be made with local ULIRGs, known to be exclusively major mergers. 
As Fig. 5 of IVeilleux et al.l (|2002l) shows, ~64% of local ULIRGs are either co alesced, or have pro jected 



nuclear separations < 2 kpc. This is what is expected from merger simulations (jMihos fc Hernquist 



19961: 



Springel et al]|2005l) . which show that a first, minor starburst is triggered upon the first encounter, and the 
most prodigious star formation rates are only achieved during the final coalescence of the two galaxies, when 
tidal torquing compresses gas into the centre and leads to a ~kpc-scale central starburst. Could it be that 
the SMG selection criterion selects a similar class of galaxies, at high redshift, and thus that the large fraction 
of compact SMGs in our sample is naturally explained as being those end-stage, coalesced major merger 
products? Indeed, as Fig.[5]shows, local ULIRGs and our SMG sample show a st rikingly similar distr ibution 
of projected separations (here we have chosen to adopt a coarser bin size than IVeilleux et al.ll2002l due to 
the comparatively small SMG sample size and the relatively large upper limits on some SMG sizes - we note 
that most of the local ULIRGs in the <6kpc bin have projected separations <2kpc). However, a direct 
comparison of the two samples is difficult. This is because Lju and Ssso^m may not always be correlated 
during a major merger, due to the way the infrared SED changes with dust temperature at these redshifts - 
for a given Ssso/im decreases wit h increasing dust tem perature. This results in a cool dust bias in the 



SMG sample, illustrated by Fig. 4 of ICasev et al.l (|2009al ). and borne out by observations (jMagnelli et al 
20101 ). This is important, since it may affect which merger stages are sampled by Lju and S^o^m selection, 



respectively: Both observations and simulations indicate that dust temperatures increase towards the final 
stages of a merger, suggesting that SMGs should sample earlier merger stages than ULIRGs. Does this imply 
that the observed ratio of binary and single SMGs is reconcilable with the single SMGs being isolated disk 
galaxies, rather than end-stage major merger products as suggested by the comparison with local ULIRGs? 
In order to answer this, we need to quantify the impact of the SMG cool dust bias on the merger stage 
sampling function. 



Veilleux et al.1 (|2002l ) use the ratio of IRAS 25 fim and 60 \xm flux as a proxy for dust temperature, and 
find that local ULIRGs with /25//60 > 0.2 ('warm ULIRGs') are only found among system with nuclear sep- 
arations below 10 kpc (no trend is evident for smaller /25//6OJ i-e. cooler systems; their Fig. 6). However, the 
percentage of systems with /25//6O > 0.2 is constant for all nuclear separations < 10 kpc. More importantly, 
only 14 out of 118 systems in their sample have /25//6O > 0.2 - i.e. less than 10% of ULIRGs are affected by 
this bias! 



Turning to simulations, Figs. 3&4 of Narayanan et al. ( 2010b[) and Fig. 1 of Narayanan et al. ( 2010a ) 



- 11 - 



show the evolution of the star formation rate, bolometric luminosity, and S^o^m during simulated gas-rich 
major mergers. These show that after the peak in star formation rate during final coalescence, Ss50/j,m appears 
to drop off more sharply than L^i, most likely due to an increasing AGN contribution effecting an increase 
in the dust temperature. However, variations between different simulations make it quite difficult to quantify 
the time a merger has both high bolometric luminosity and low Ss50fj, m (D. Narayanan, priv. comm.). 

Another way to tackle this question may be to compare the AGN luminosity contributions in local 
ULIRGs and SMGs, since from t he above we would e xpect SMGs to be biased against the merger phases 
with strong AGN contributions. IVeilleux et al.l ([20091 ) find local ULIRG luminositie s to have ~ 39 % AGN 



contribution; estimates for SMGs are more uncertain, but span a range of 8-35% (Swinbank ct al. 2004; 



Alexander et al.ll2005[Eo08aHLutz et al.ll200d : lvaliante et al.ll2007tlPope et al.l2008l ; [Menendez-Delmestre et al 
2009). This seems to confirm a significant bias - however, one must bear in mind the differences between 



high-z and local galaxies; their absolute AGN luminosities should be comparable (similar black hole masses), 
but the stellar lum i nosities are factor s of several larger in high-z galaxies (due to the larger gas fractions, 
Tacconi et al.ll201Clt iDaddi et al.ll2010l ). and we therefore expect significantly smaller relative AGN contribu- 
tions in high-z galaxies. Lower relative AGN contributions in SMGs compared to local ULIRGs are therefore 
not necessarily a result of a merger stage bias. 

The cool dust bias of the S'sso/im criterion implies that there should be a population o f hot du st, ultra- 
/hyperluminous high-z galaxies which are not detected as SMGs. And indeed, I Casey et alj (|2009af) through 
70 /im selection find a population of z > 1 galaxies with an average dust temperature of 52 K, which are not 
detected at 850 //to. These systems have a mean far-infrared luminosity of 1.9 x lO 12 ^©, only slightly lower 
than that of SMGs, and thus it is plausible that they consitute the hot dust extension of the SMG sample. 
Although uncertain, the number density of these galaxies is estimated to be approximately five times lower 
than that of SMGs - implying that, if indeed all these systems are major mergers, a high-z major merger 
may be undetectable as an SMG for approximately l/6th of its lifetime; or «25% of the time spent in the 
'final coalescence' phase, if all these hot dust systems are in that stage (as one would expect, based on the 
arguments above). This is larger than, but consistent with, the ~ 10% of local ULIRGs showing a dust 
temperature - merger stage correlation. 

It is clear that the merger stage sampling bias of Ssso^m vs. Lib. selection is difficult to quantify. We 
can only say with confidence that Sgso^m sampling of a major merger will likely result in a somewhat smaller 
fraction of coalesced systems, but it certainly will not exclude all of them. Loca l major mergers, sel ected 
via Ljb, consist to ~64% of systems with projected nuclear separations <2kpc ( Veilleux et al. 2002); our 
sample of 12 SMGs contains 7 sources which are either poorly resolved or clearly single galaxies (Table [2]). If 
we conservatively assume that a major merger is undetectable as an SMG for 50% (30%) of the time spent 
in the 'final coalescence' stage, this implies we would expect the fraction of coalesced systems in the sample 
to decrease from ~ 64% to about 32% (45%). Thus at most 3(1) out of the 7 compact/coalesced sources in 
our sample could plausibly be isolated disk galaxies, rather than end-stage major mergers. 



This corroborates the argument made by iTacconi et al.l (|2006f ) , that the extreme surface densities of 
the compact SMGs can plausibly only be achieved during a major merger, an d is further supported by our 
results of ij4.2.2[ where we show that our confirmation of the iHainlind ([20081 ) stellar mass estimates agrees 
with cosmological models in which SMGs result from major mergers. 



These three independent arguments all point in one direction, namely that most, if not all, SMGs are 
major mergers. 
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5. Summary and Conclusions 



Our investigation of twelve sub-arcsecond resolution interferometric CO line observations of SMGs lends 
strong support to a merger origin of a large majority of this class of galaxies. 

• Five of our systems consist of two spatially distinct galaxies, and the remaining seven have either 
disturbed morphologies typical for advanced, pre-coalescence mergers or are compact, dense galaxies 
which plausibly are late-stage, coalesced mergers. 

• All four systems consisting of two resolved galaxies in which the masses of the two separate components 
could be derived, have mass ratios of 1:3 or closer, within the errors. This demonstrates that they are 
major mergers. 

• Our gas and dynamical mass estim ates are physica lly consistent w ith the lower end of the range of stellar 



masses in the literature, those of Hainline ( 20081) adjusted to a Chabrier ( 2003 ) IMF. This supports 



hierarchical merger models, whereas it disfavours cold accretion flow models of SMG formation. 

• The majority (~64%) of local infrared-luminosity-selected major mergers (ULIRGs) are coalesced or 
have nuclear separations < 2 kpc. This suggests that the large fraction of binary SMGs in our sample 
also implies a major merger origin for at least some of the compact/coalesced SMGs. Due to the 
uncertain difference in merger stage sampling of Lm > 10 12 L© and Sssofim > 5 mJy selection, we 
can only estimate the corresponding ratio of binary/coalesced mergers expected in our SMG sample. 
However, even very conservative assumptions indicate that at most 3 of the 7 compact/coalesced 
sources in our sample could plausibly be isolated disk galaxies, rather than end-stage major mergers. 

• SMG size measurements (Gaussian FWHM) from radio and CO fluxes are both consistent with diam- 



eters ~ 5 kpc (somewhat larger than the results of iTacconi et a 



number statist i cs). T hese sizes are in agreement with work by 



( 200 61) - most plausibly due to small 
" (|2009l ) and 



Menendez-Delmestre et al 



Hainline et al.l (]2009f ). who find that the star formation in SMGs is more extended than the kpc-scale 



central starburst in local ULIRGs. This is consistent with what is expected from theoretical pred ictions 
and simulations of mergers with higher gas fractions (e.g. Mihosl 1999 : Narayanan et al. 2010b[ ). 



Theory and observations are only beginning to uncover how differently galaxy evolution proceeds at 
high redshifts, with the 'cold accretion flow' framework successfully addressing many of the shortcomings of 
the hierarchical merging scenario. Particularly in the lower-luminosity tail of the high-z galaxy population, 
future observational and theoretical efforts are clearly necessary to fully clarify the picture. However, the 
evidence presented here overwhelmingly supports a framework in which the most luminous z ~ 2 galaxies are 
powered by major mergers. 



We would like to thank the anonymous referee for a thoughtful and thorough reading of the manuscript, 
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discussion of his simulations and the merger evolution of Ssso^m vs. Lir. We thank the staff of the IRAM 
observatory for their support of this program. The Dark Cosmology Centre is funded by DNRF. TRG 
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Table 1. Measured & Derived Quantities 



salaxv 


obs. CO 


Sco^V 


HWHM UV 


L W Jl IVl 


M gas 


J Q,yn 




transition 


[Jykms -1 ] 


fkpcl 

L I J 


[kms~ 1 


\10 10 Mp.} 


[10 10 M CT 1 


SMMJ09431, SW 


6-5 1 


1.8±0.3 


1.13±0.23 


390±45 


7.10±2.44 


4.81±1.26 




4-3 3 


1.1±0.2 


1.85±0.60 


400±45 


5.64±1.98 


8.29±2.99 




1-0 6 


<0.13 


n/a 


n/a 


<8.67 


n/a 


SMM J09431, NE 


6-5 1 


1.1±0.2 


n/a 


n/a 


4.34±1.53 


n/a 




4-3 3 


<0.4 


n/a 


n/a 


<2.05 


n/a 


SMM J131201 


6-5 1 


2.54±0.1 


4.09±0.51 


911±140 


13.43±4.97 


95.04±23.81 




4-3 5 


1.7±0.3 


n/a 


530±50 


11.67±4.06 


n/a 




1-0 6 


0.42±0.07 


n/a 


1040±190 


28.83±4.81 


n/a 


SMMJ105141, SW 


4-3 1 


1.53±0.11 


<1.40 


436±50 


1.63±0.51 


<7.45 


SMM J105141, NE 


4-3 1 


5.87±0.19 


1.77±0.14 


383±50 


6.72±1.89 


7.27±1.46 


HDF169 


4-3 1 


2.25±0.11 


2.15±0.24 


431±60 


2.45±0.75 


11.18±2.53 




2-1 8 


3.45±0.93 


<16.5 


560±90 


12.47±5.03 


<145 


HDF76 


6 _ 5 2,3 


2.3±0.4 


0.95±0.4 


600±80 


5.76±2.00 


9.58±4.42 




3-2 2 < 3 


1.6±0.2 


<2.1 


600±50 


7.44±2.42 


<21.17 


HDF242, SW 


3-2 1 * 2 ' 3 


0.59±0.08 


2.51±0.86 


426±60 


3.41±1.13 


12.75±5.10 


HDF 242, NE 


3-2 1 ' 2 ' 3 


0.32±0.07 


3.45±0.98 


394±60 


1.85±0.68 


15.00±5.35 


N2850.2 


7 _ 6 2,3 


3.3±0.5 


0.8±0.5 


790±50 


n/a 


13.98±8.83 




3-2 3 - 5 


1.8±0.2 


n/a 


800±50 


10.11±3.23 


n/a 


N2850.4, total 


3-2 2 < 3 


2.3±0.3 


2.48±1.24 


710±50 


12.28±4.02 


35.00±17.85 


N2850.4, SW 


7 . 6 1,2,3 


2.01±0.18 


0.54±0.37 


373±50 


n/a 


2.10±1.50 


N2850.4, NE 


7 . 6 1,2,3 


1.36±0.18 


<1.70 


139±50 


n/a 


<0.92 


SMM J044307 


3-2 3 


1.4±0.2 


<3.26 


350±40 


8.20±2.73 


<11.18 


SMM J02399 


1-0 4 


n/a 


n/a 


n/a 


100±30 


n/a 


HDF 132 


4-3 7 


1.27±0.14 


2.65 


n/a 


3.49±1.12 


21.80±1.80 


Lockman 38 


3-2 7 


0.63±0.08 


4.75 


n/a 


1.49±0.49 


6.20±1.10 



1 this work 
Tacconi et al. (2008) 



Tacconi et all l|2006l) 



z Ivison et all (l2010h 
Greve et all d2005h 



lHainline et al.l (|2006f ) 
Bothwell et all <|2010| ) 
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Fraver et al.l ([2008I ) 



Note. — Measured quantities for the twelve systems in our sample. New observations presented in this 
paper are in bold face. If a system consists of two galaxies, SW or NE indicates which galaxy is referred 
to. As outlined in the text, the HWHM of Gaussian fits to the source in the iro-plane were taken as 
best esti mates of the intrinsic R m?., except for HDF132 and Lockman38, for which the half light radii as 
quoted bv lBothwell et al.l (|2010h are given. The quoted values for dynamical masses were calc ulated using the 
' isotro pic virial estimator' (see i)3.1l for details), except for HDF 132 and Lockman38, for which iBothwell et al. 
(|2010l ) estimated the quoted masses using disk rotatio n mode l s. Ga s masses were calculated by converting 
the measured line flux to CO(l-O) flux using Fig. 3 of IWeissI (|2007l) . assuming a ULIRG conversion factor 
a = lM /Kkms~ 1 pc 2 , and including a factor 1.36 correction for hel ium. For SMM J09 431, measurements 
were corrected for gravitational lensing magnification of factor 1.3 (|Cowie et al.l 12002). If a source was 
unresolved, the beam HWHM was taken as an upper limit to R\/2- 
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Table 2. Merger Characteristics 



galaxy 


z 


proj. sep. 
[kpc] 


mass ratio (gas) 


mass ratio (dyn) 


SMMJ09431 1 - 3 


3.35 


30±3 


0.41±0.2 


n/a 


SMMJ131201 1 


3.41 


<6 


n/a 


n/a 


SMMJ105141 1 


1.21 


4.9±0.5 


0.26±0.02 


>0.98 


HDF 169 1 


1.22 


<1.5 


coal. 


coal. 


HDF 76 2 ' 3 


2.20 


<2 


coal. 


coal. 


HDF 242 1 ' 2 - 3 


2.49 


22±2 


0.54±0.14 


0.85±0.49 


N2850.2 2 ' 3 


2.45 


<3 


coal. 


coal. 


N2850.4 1 ' 2 - 3 


2.38 


3.7±0.4 


0.68±0.11 


<0.44 


SMMJ044307 3 


2.51 


<4 


coal. 


coal. 


SMMJ02399 4 


2.80 


-8 


n/a 


n/a 


HDF 132 5 


2.00 


<3 


coal. 


coal. 


Lockman 38 5 


1.52 


<4 


coal. 


coal. 



Hhis work 
Tacconi et all |2008l ) 



foacconi et all |2006l ) 
1 Ivison et all (|2010l ) 



Bothwell et all (|2010h 



Note. — Projected separations and mass ratios. Coal, denotes coalesced 
system. The distribution of projected separations is also plotted in Fig. [5] For 
SMM J09431, the CO(6-5) data presented in this work were used to derive the 
gas mass ratio. For SMM J02399, no information on dynamical or gas masses 
of the subcomponents was available in the literature. As can be seen, the 
mass ratios are all consistent with major (mass ratio 1:3 or larger) mergers, 
within the error ranges. 



Table 3. Mass Measurements 



galaxy 


M gas 


Mdyn 


M+ 1 


M* 2 


Mj 




[1O 1O M ] 


[1O 1O M ] 


[1O 1O M ] 


[1O 1O M ] 


[1O 1O M ] 



SMMJ131201 


18.0±9.4 


95.0±23.8 


n/a 


10.2lt 9 3 


13.5 


HDF 169 


7.5±5.0 


11.2±2.5 


n/a 


3 2+ 11 


26.9 


HDF 76 


6.6±1.6 


9.6±4.4 


12^ 




28.2 


HDF 242 (SW + NE) 


5.3±1.3 


27.8±7.4 


i 9+2.5 
lz -4.5 




52.4 


N2850.2 


10.1±3.2 


14.0±8.8 


25±7.5 


q 9 + 1.3 
°- z -0.9 


34.6 


N2850.4 (SW + NE) 


12.3±4.0 


35.0±17.9 


231^ 


17.4™ 


21.9 


HDF 132 


3.5±1.1 


21.8±1.8 


n/a 


2-oo±8;g 


29.5 


SMMJ16371+4053 


5.3±1.1 


n/a 


n/a 




47.8 


SMMJ22174+0015 


6.7±1.7 


n/a 


n/a 


3-4±?i 


13.8 



Tacconi et al.l (|2008l ) 



feainlinel |2008l ) 
^Michalowski et al.l |2009l) 



Note. — Stellar, gas, and dynamical mass me asurements for the seven galax ies which 
are co ntained both in our sample and those of Michalowski et al. ( 20091) and Hainline! 
( 20081 ) , and two additional galaxies contai ned in the latter tw o samples for which we cal - 
culate gas masses using fluxes measured by Greve et al. ( 20051) . Michalowski et al. ( 2009 ) 
do not quote un certainties for their stellar masses. Stellar masses were corrected for a 
Chabrierl (|2003l ) IMF. Binaries are treated as one system. Where multiple transitions 



were observed, dynamical and gas masses are averages of those deriv ed for the different 



transi tions, see Table[T] As can be seen, the stellar m a sses d erived bv lMichalowski et al 



( 20091) are significantly larger than those of iHainlinej (|2008l ). with only the latter being 
in good agreement with our gas and dynamical mass measurements. 
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Fig. 1. — SMMJ09431; this system consists of two galaxies separated by ~30kpc. The south-west one 
(H7) is a compact, rotating, dispersion-dominated disk galaxy. The left spectrum is of the north-east galaxy 
(H6), the right spectrum of the south-west galaxy. On the right are, top to bottom, south-west galaxy flux, 
velocity, and dispersion map, on the left north-east galaxy flux (the difference to the large-scale map being 
due to the H6 map having been integrated over the observed line extent rather than the full wavelength 
range) and integrated spectrum. The beam size (0.64"x0.48", P.A. = 13.5°) is displayed in the lower left 
corner of the flux maps. North is up and east is to the left. Scale bars denote 1". 



SMMJ131201 
, CO 6-5 





Fig. 2. — SMM J131201; this system displays a complex, irregular morphology, it most likely is an advanced 
pre-coalescence merger. The RGB-image (mapping red, green, and blue parts of the spectrum) shows no 
indication of rotation. Spectra are, left: north-east (top) and south-east (bottom) arm, right: entire system 
(top), and central part (bottom). The beam size (0.59"x0.47", P.A. = 50.9°) is displayed in the lower left 
corner of the flux maps. North is up and east is to the left. Scale bar denotes 1". 
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Fig. 3. — SMM J105141; as the RGB-map (top; mapping red, green, and blue parts of the spectrum) clearly 
shows, this system consists of two galaxies that are separated in velocity by ~ 400 km s _1 . Top right spectrum 
is of 'main' component (north-east galaxy), bottom right spectrum of the blue-shifted component (south-west 
galaxy). As can be seen from the velocity and dispersion maps (left, middle and bottom), the north-east 
galaxy is a dispersion-dominated disk galaxy. The beam size (0.46"x0.34", P.A. = 28.7°) is displayed in the 
lower left corner of the flux maps. North is up and east is to the left. Scale bar denotes 1". 
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Fig. 4. — HDF 169; the position-velocity diagram (taken along the dotted line shown in the flux map, 
from lower right to upper left) and RGB-map (mapping red, green, and blue parts of the spectrum) indicate 
rotation, which is also supported by the integrated spectrum which appears double-peaked. This galaxy thus 
could either be a disk galaxy or a compact, rotating coalesced merger remnant. The beam size (0.48" x 0.35", 
P. A. = 161°) is displayed in the lower left corner of the flux maps. North is up and east is to the left. Scale 
bar denotes 1". 
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Fig. 5. — Dist r ibutio n of projected separations. Red: SMGs, this work, 12 systems. Blue: local ULIRGs, 



Veilleux et al.l (|2002l ). 118 systems. Note that the last bin contains all systems with projected separations 



larger than 24 kpc. We note that for the SMGs, of the 9 systems in the first bin (projected separation 
smaller than 6 kpc), 7 are either poorly resolved or coalesced, and 2 are clearly recognisable as close binaries; 
the SMG sample thus divides into 5 well-resolved binary systems and 7 which are either poorly resolved or 
clearly compact, rotating systems. It should be noted that ~ 64% of local ULIRGs have projected separations 
< 2 kpc, however due to the larger upper limits of some of the SMGs, and the increased uncertainties of a 
finer resolution (due to our small sample size) , we have opted to use a coarser resolution. It is evident that 
SMGs and local ULIRGs display a very similar distribution of projected separations. 
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Fig. 6. — Plot of (M gas + M*) vs Md vn for the seven gal axies for which all three mass estim ates are availabl e 
(Table [3]). Blue cross es denote 
both corrected for a 

stellar masses. The dashed line demarcates M a 



Mich alowski et al. (2009) st ellar masses , red squares those of lHainlind ((2008) ; 



Chabrierl (|2003l ) IMF. I Michalowski et al.l (|2009l ) do not quote uncertainties for their 

M* = Mdyn] in order for the mass estimates to be 



physically realistic, galaxies should lie on or below this l i ne. W ithin the errors, the iHamlme] |2008l ) stellar 
ma sses pass t his te st, whereas those of iMichalowski et al. I (|2009h do not. This constitutes strong support for 
the lHainlinel (|2008() results. 



